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Abstract 



Wc consider the diffusion approximation of branching processes in random environment 
(BPREs). This diffusion approximation is similar to and mathematically more tractable 
than BPREs. We obtain the exact asymptotic behavior of the survival probability. As in the 
case of BPREs, there is a phase transition in the subcritical regime due to different survival 
opportunities. In addition, we characterize the process conditioned to never go extinct and 
T^ • establish a backbone construction. In the strongly subcritical regime, mean offspring numbers 

C^ , are increased but still subcritical in the process conditioned to never go extinct. Here survival 

is solely due to an immortal individual, whose offspring are the ancestors of additional families. 
In the weakly subcritical regime, the mean offspring number is supercritical in the process 
conditioned to never go extinct. Thus this process survives with positive probability even if 

►^ . there was no immortal individual. 

m 

t — , 

l^ . 1 Introduction and main results 

(N 

t^^ . Branching processes in random environment (BPREs) have been introduced by Smith and 

^^ i Wilkinson (1969) (see also Smith (1968)) and have attracted considerable interest in the last 

' decade (e.g. [aHSl El [El [2011211 123 SHI HI])- On the one hand this is due to the more reahstic 

model compared with classical branching processes. On the other hand this is due to interesting 

properties such as phase transitions in the subcritical regime. Here we consider the diffusion 

K> ' approximation of BPREs which can be viewed as continuous mass branching process in ran- 

j^ ' dom environment. Our resuhs are quahtatively analogous to discrete mass BPREs. The main 

observation of this article is that the diffusion approximation of BPREs is a simple model (3 

parameters) having explicit formulas for various expressions. In particular, the contributions of 

the branching process and of the environment are explicit in terms of the parameters. These 

properties make the diffusion approximation interesting for applications. 
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The diffusion approximation of BPREs has been conjectured by Keiding (1975) and has been 
estabhshed by Kurtz (1978). This diffusion approximation is the strong solution {Zt,St)t>o of 
the stochastic differential equations (SDEs) 

la^.Ztdt + ZtdSt + J^tdW^^^ 
t = adt+ y^dW}' 



for t > where 5o = 0. The parameters satisfy a £ ]R,,ae € [0,oo) and df, G (0, oo). The 
processes [WJ^ )t>o and (W^ )t>o are independent standard Brownian motions. We denote the 
process {Zt,St)t>o as branching diffusion in random environment (BDRE). Moreover, we will 
refer to {St)t>Q as the associated Brownian motion, which is a non-standard Brownian motion. 
To be accurate, the conjecture of Keiding (1975) did not include the term ^a'^Ztdt which is 
a characteristic part for random environment as we will see below. Moreover, Helland (1981) 
gave an inaccurate "proof" of N. Keiding's conjecture. So we attribute the correct statement of 
the diffusion approximation of BPREs to Kurtz (1978). The BDRE has not been studied since 
Kurtz (1978). For this reason, we first discuss properties of the BDRE obtained by Kurtz (1978) 
beginning with the diffusion approximation. 

First we introduce BPREs in order to state the diffusion approximation. Our formulation 
follows the notation of Afanasyev et al. [1]. Let A be the Polish space of probability measures on 
No := {0, 1,2,.. .} equipped with the metric of total variation. Fix n G N := {1, 2, . . .} for the 
moment. Let II^"' = {Qi ,Q2 ■, • • •) be a sequence of independent and identically distributed 
random variables taking values in A. Conditioned on II'") the BPRE [Z^^^ )iQ-^^ is defined 
recursively through 

4:1 -E^?' ^^No, (2) 

i=i 

where Zq is independent of H^""' and where (Cii )i,jeNo conditioned on n^") are independent 
random variables with distribution 

P(eg) = A;|nW) = qJ")(A:), yj,i,k G No. (3) 

Let the mean of the environment at time i G N be defined through 

oo 

m{Q\-^):=Y,kQ[-\k). (4) 

fc=0 

Define a continuous time version of the BPRE through Z^ := Z,", where [t\ := max{m G 
No: rn < t} for every t G [0, oo). The associated random walk (5j )^>q is defined through 

Lij-i 
^J") := ^ E l°g ("^(QI"^))' i e [0,oo). (5) 



This random walk is central for the BPRE as it determines the mean of the BPRE: 



W-i 



E[Z,(")|n(")] =E[zi")] H m(Qf^) =E[zi")]exp 



S, 



(n) 



i=l 



t G [0, oo) 



(6) 



We included the factor ^/n in the definition of the associated random walk to have the usual 
scaling in the limit as n — ?• oo. 

Next we let n — t- oo to obtain the diffusion approximation. The following assumptions mainly 
ensure that the associated random walk converges to a Brownian motion with infinitesimal drift 
a S R and infinitesimal standard deviation ae G [0,oo): 



lim n • E 

n— >oo 



m(Ql"))-l 



lim n ■ E 



sup E 

nSN 



fc=0 



m 



(Qi' 



C«)^ 



{m{Q\ 

3 
- 1 



(")^ 



ale [0,oo) 



i(") 



QT'ik) 



< oo. 



(7) 
(8) 

(9) 



So the branching process is near-critical as mi^Q^^ ) — )• 1 in distribution as n — )■ oo. If o"e > 0, then 
the environment comprises both supercritical and subcritical phases. In addition, we suppose 
that the mean branching variance in each generation stabilizes in the limit as n — t- oo 



limE 

n— >oo 



oo 

[E 

fc=0 



m 



^ lfQ^l^\k)]=aU{Q,^). 



{Ql 



(10) 



So a is a parameter of expected super-/subcriticality, ae is a parameter for the standard deviation 
of the offspring mean around the critical value 1 and a"^ is the mean offspring variance per 
individual per generation. Under the above assumptions, Kurtz (1978) has shown in his Corollary 
2.18 that the suitably rescaled BPRE converges in distribution to a diffusion. 

Proposition 1. Assume that Z^ jn — >• Zq in distribution as n — >• oo. Under the assump- 
tions dZ}, dSD, Q and ([TOD we have that 






n 



(Zi,5i)^>Q 



:ii) 



in the Skorohod topology (see e.g. \17^ ) where the limiting diffusion is the strong solution of the 
SDEs ^. 

The assumptions of Corollary 2.18 of [3U] are checked in Section [21 

Inserting the random environment {St)t>o into the diffusion equation ([1]) of the BDRE, we 
see that {Zt)t>o solves the stochastic differential equation 



dZt =(a + ]^al^ Zt dt + ^ alZl dW^"-^ + .plZtdwf^ 



(12) 



for t € [0,00). Comparing with Feller's branching diffusion (i.e. ()12p with o"e = 0) there are 
two differences. First there is an additional drift term ^a'^Ztdt. Second there is an additional 

diffusion term CeZtdWf . Both terms originate in the conditional expectation 

E[Zt\St]=E[Zo]exp{St), tG[0,oo), (13) 

almost surely, which is a geometric Brownian motion and solves the SDE 

dYt = (a + ^^2) Ytdt + aeYt dw}'^ , Yo = E [Zq] , (14) 

for t > 0. 

Now we come to properties of the BDRE which embody the branching property conditioned 
on the environment. Theorem 2.10 of Kurtz (1978) implies that the BDRE is in fact a reweighted 
and time-changed branching diffusion (see Lemma [TBI below for a different proof): 

Proposition 2. Assume a G R, at G (0, 00) and a^ G [0,oo). Let (W^ )t>o and (VFj )t>o be 
independent standard Brownian motions. Let (i*t)t>o be the strong solution of 



dFt = ^/FtdW^^^ (15) 



r{e) 



fort£ [0,00) and let St := at + agWf fortG [0,oo). Moreover define the random time {T{t))t>o 
through 

T{t) := / e-'^'a^ds (16) 



for t G [0, 00). Then 

iFrit)e'\St)^^^ (17) 

is a weak solution of ^, that is, is a version of the BDRE ([T]). 

Due to this property, many results on Feller's branching diffusion carry over to the BDRE ([T]). 
For example, (Z^e" *)j>o is a time-changed Feller branching diffusion and is therefore infinitely 
divisible. Another simple implication of Proposition [2] is an explicit formula for the Laplace 
transform of the BDRE ([1]) conditioned on the environment, which has not been reported yet. 
We agree on the convention that 

c I 00 if c G (0, ool c „ r s , . -. 

- := { , — := for c G 0, 00) and that • 00 = 0. (18) 

lO ifc = 00 

Throughout the paper, the notation P^ and E^ refers to the starting point of the involved process, 
e.g. F'{Zt G •) := F{Zt G -IZq = z) for z G [0, 00) or p(^'^)((Zi, 5*) G •) := P((Zi,5i) G 
•|(Zo,S'o) = iz,s)) for z G [0,00) and s G R. 

Corollary 3. Assume a G R, Cf, G (0, 00) and Ue G [0, 00). Let {Zt,St)t>o be the strong solution 
of ^. Then we have that 

E»'' h (-...) I (5.). J . exp ( - ,^^^,_',^ ) (10) 



/g ^ exp [- Ss)ds + J exp(-5t) - 



for all t, z, X £ [0, 00) almost surely. 



The proof is deferred to Section [3l If a^ = 0, then (|19p is just the Laplace transform of Feller's 
branching diffusion with criticality parameter a and branching rate <t^. 

The simplicity of the right-hand side of ()19p derives from the fact that the distribution of the 
integral of the geometric Brownian motion with drift /? G IR, 

^f^= / exp(2(/3s + W^i^)))ds, tG[0,oo), (20) 

Jo 

is well understood (see e.g. |10 p i4 1[33ll42| ). Even more, the density of the joint distribution of 

{A^ ' , W^ + (3t) is known rather explicitly for every t G [0, oo). Define 

at{x, u) du := P(Af ^ G du\W^^^ + /3t = x) (21) 

for t,ue (0, cx)) and x G R. Then the density of {Af\ WJ;''^ + f3t) satisfies that 
1 / X^N , , 1 /I 



V27ri 
where 



exp [-^)at{x, u) = -exp (^- — (1 + e^-)]^,./^*) (22) 



^''^^^ ^ /K-r ^^P V 2t / / '^^^ '-" 2t '^ ^^^ '-"^ cosh(y)) sinh(y) sin (^— j dy (23) 



for all t,u,r G (0, oo) and x G R, see Proposition 2 of Yor (1992). Using the explicit formula (|22p 
allows to answer rather fine questions by elementary (but sometimes nontrivial) calculations. 
Thereby the BDRE ([TJ becomes one of the most tractable processes in the class of BPREs. 

The following corollary of Corollary [3] provides an explicit expression for the survival proba- 
bility. Define the parameter /? G [—00,00] and the function /: [0,oo] — )■ [0, 1] through 

/3 := 2 ^^d f{x) -=1 - exp( 1 • a; ) , xG[0, 00]. (24) 






Corollary 4. Assume a G R, Cf, G (0, 00) and ae G [0, 00). Let {Zt,St)t>o be the strong solution 
of (dD. Then 

F(--°)(z,>0|(^.).,,)=l-exp(- - ) (25) 

^ io^exp(-5,)ds/ 

for every t G (0, 00) and every z G [0,oo) almost surely. If j5 > —1 and ae > 0, then 

P^(Zi>0)=E[/(^^)] = r f{za)pt„2j^^p{a)da (26) 



^9 4^ 







for every t G (0, 00) and every z G [0,oo) where the density function ofl/{2A\, ) satisfies 
Py^i3{a)da := F\^ ^ G daj 



00 roc 






r2 

sinh(0 cosh(,^) sin(7r^/t') 
oT^ dt, ds da 

(s-f(cOsh(0)2) 2 



T[P±iy-a-^P+^^'^ / e-«'/2-^(/3-i)/2g-a. (27) 



on (0, cxo) for every v € (0, cxo). 



The proof is deferred to Section [31 

The asymptotic behavior of the survival probabihty strongly depends on a. As in the case of 
classical branching processes, the survival probability stays positive, converges to zero polynomi- 
ally fast or converges to zero exponentially fast according to whether the process is supercritical 
(a > 0), critical {a = 0) or subcritical {a < 0), respectively. Now in case of a random environ- 
ment it is known for BPREs that there is another phase transition in the subcritical regime. For 
BDREs this phase transition turns out to occur at a = — CTg. We adopt the standard notation of 
the literature on BPREs for the different regimes and say that the BDRE is weakly subcritical 
if — (jg < a < 0, intermediately subcritical if —a'^ = a and strongly subcritical if a < — o"g. 
The following theorem establishes the asymptotic behavior of the survival probability of BDREs 
including explicit expressions for the limiting constants. 

Theorem 5. Assume a G R and ah^cTe G (0,cx3). Let {Zt,St)t>o be the strong solution of ([1]). 
Then we have that 



hm ¥^Zt>0 

t—>-oo 



limVtF^ Zt >0 



^3 ^t 

lim vi e^'^e 



lim V i e 2 

t—^oo 



t mz 



lim e 

i— >oo 



-("+4) 



Zt>0 
^Zt>0 

t,az( ^^ >0 



2a 

l + ^-z 



>0 



V2 



log(l + ^- 
V at 



z >0 



-^ / f{za)(l)fs{a)da > 

^i Jo 

Z ^ n > 



y/na_ 



z2^^^>0 



if a> 


(28) 


ifa = 


(29) 


^/4G(-l,0) 


(30) 




(31) 


^/"<-i 


(32) 



for every z G (0, oo), where /3 := — % and where (pf^: (0,oo 
for every a G (0, oo). 



r(^rf)e-'^a-^/2^('^-i)/2g-« 



— )■ (0, c«) is defined as 

sinh(^) cosh(^)^ 

j+^djau 

{u + a(cosh(^))2) 2 



(33) 



The proof is deferred to Section HI 

Let us compare the convergence rate of the survival probability with the classical case o"e = 
of Feller's branching diffusion. Pars pro toto we discuss the critical regime a = 0. In that 
case, the survival probability of Feller's branching diffusion is of order 0(j) whereas it is of 
order O(^) if Ue > as t — )• oo. So a branching process in random environment has a higher 
probability to survive. The reason for this is that there is a positive probability of experiencing 
a long supercritical phase. More precisely, for every e > 0, the event that the critical Brownian 
motion {Ss)s>o stays above e from time e until time t is of order O(^) as t — t- oo. On this event 



the branching process is supercritical and survives with positive probabiUty. This explains the 
slower convergence order O(^) ast— t-oo. 

2 

Note that the expectation W[Zt] = zexp ((a + ^)t), z S (0, oo), changes its qualitative 

2 

behavior asi— )-ooata = — ^. The phase transition for the survival probability, however, is at 
a = —(j1. Here is an heuristic. If the associated Brownian motion (drift a < 0) is negative for 
almost all of the time, then we expect the BDRE to behave like Feller's branching diffusion. In 

2 

that case we expect that P^(Zt > 0) ~ const ■ W{Zt) = const ■ zexp ((q + ^)t), z G (0, oo), as 
t — )• oo. This gives indeed the exponential decay rate in the strongly subcritical regime. However, 
the associated Brownian motion might be positive until time t > 0. On this event the BDRE is 
supercritical and survives with positive probability. The probability of this event decreases like 
exp (— ^T^) (times polynomial terms) as t — )• oo. This exponential decay rate follows from an 
application of the Cameron-Martin-Girsanov theorem (e.g. Theorem IV. 38. 5 in [35])- This gives 
the exponential decay rate in the weakly subcritical regime. Now the phase transition for the 

2 2 

survival probability occurs when these two exponential decay rates a + ^ and — ^ coincide, 

namely at a = — o"g. 

For BPREs Afanasyev (1979) was the first to observe different regimes for the survival proba- 
bility in the subcritical regime. Independently hereof Dekking (1987) rediscovered this dichotomy. 
For more recent results on the speed of decay of the survival probability, see Corollary 1.2 of [3] 
for the critical case. Corollary 1.2 of [3] for the weakly subcritical case. Theorem 1 of [IQ] for 
the intermediately subcritical case and Theorem 1.1 of [5] for the strongly subcritical case). Its 
derivation, however, is sometimes involved and, in general, there are no simple expressions for the 
limiting constants. Only the case of linear-fractional offspring distributions is known to admit 
explicit limiting constants, see |2]. 

Next we investigate the event of survival in more detail and condition the BDRE on the event 
of ultimate survival. The method of conditioning a Markov process to stay nonnegative has been 
applied in various situations (e.g. [Z1E1131]). However, we have not found a suitable formulation 
for the case of multi-dimensional diffusions. As such a formulation is of independent interest, we 
include it in the following lemma. For this, define a set Q of functions as 

Q := \q: [0, oo) -^ [0,oo) lim ^i^i-fl = 1 for all s G [0,oo) 1 . (34) 

Note that g G Q if and only if g o log is slowly varying at infinity (see Galambos and Seneta 1973 
for this notion). 

Lemma 6. Let d,m ^ ¥i and let I C R and A C I be Borel measurable sets. In addition let 
^: / — 7- R and a: I ^ R"^^™ be Borel measurable functions. Moreover, let {Xt)t>o be a Markov 
process and a weak solution of the stochastic differential equation 

dXt = ii{Xt) dt + a{Xt) dWt (35) 

for t G [0, oo) with initial value Xq G / where {Wt)t>o is an m-dimensional standard Brownian 
m,otion. Assume that there exist a twice continuously dijferentiable function r]: I ^ [0,oo), a 
function q £ Q and values X,p £ [0,oo) with the following properties: 



lim- 



■t~>-oo ' 



q(^iy\t p^ Xt ^ = r/(x) for all xe I, 



• E^[||Xip] <oo for allxel and t e [0,oo). 
Define I := {x £ I : r]{x) > 0}. Then there exists a process {Xt)t>o with state space I such that 



limP^ 

T-s>oo 



Xt^A 



.^*)s6[0,t] ^ 



(^s)se[o,t] 
for all X £ I and all t G [0,oo), such that {Xt)t>o is a weak solution of the SDE 



dXt 



aa 



{Xt) dt + fiiXt) dt + a{Xt) dWt 



for t G [0, oo) where V = I -^, . . . , -^ I and such that 



E^ HXt)] 



E- [viXt)g{Xt)] 



-xt 



r]{x) 



(36) 



(37) 



(38) 



for all X £ I, t £ [0, oo) and all Borel measurable functions g: I ^ [0, oo). 

The proof is deferred to Section [SI 

We will apply Lemma [6] to om' bivariate process {Zt,St)t>o and to the set A = {0} x R. 
LemmaElshows that the limiting constants of Theorem [S] play an important role for conditioning 
on ultimate survival. Theorem [5] shows that r]{z,s) = 'i?(z), {z,s) £ [0, oo) x R, where the 
function t9: [0, oo) — t- [0,oo) is defined through 

2a 

1 - (l + fl • ^) "^ if a > 
^logfl + ( 



^(z) 



-^ log 1 + % • z if a = 
j,!^f{za)<Pp{a)da if ^g (-1,0) 



(39) 



^% > if 4 = -1 



y/lTO 



z2 



\i^<-l 



for every z £ [0, oo) where /3 and (pp are defined in Theorem [5l Moreover, define A := for a > 0, 

2 2 

A := Y^ for a £ (— UgjO) and A := — (a+ ^) for a < — 0"g. The following theorem characterizes 
the BDRE conditioned to never go extinct. 

Theorem 7. Assume a G IR and cri,,ae £ (0,oo). Let {Zt,St)t>o be the strong solution of ([T]). 
Let {Zt,St)t>o denote the process {Zt,St)t>o conditioned to never go extinct. Then this process 
is a weak solution of the SDE 



dZt 
dSt 



ai^dSzt + lalZt ] dt + ZtdSt + ./^dl^f ^ 



^Zt 



a + alVS-Zt] dt + aedwi'^ 



(40) 



HZt 



for t G [0, cxo) and satisfies 

e(^'^) [g {Zt, St)] = ^lE^^'^) WiZMZt, St)] (41) 

for all (z, s) G (0, oo) x R, t G [0, oo) and all Borel measurable functions g : [0, oo) x K, — )• [0, oo). 
If a £ (— o"g,oo), then the function (0, oo) B z ^ alzd {z)/'d{z) G IR, is strictly monotonic 
decreasing and satisfies YwUz^q a^z-d {z)/'d{z) = a1 and lim^_>.oo alz-d {z)/'d{z) = max(— a,0). // 
a G (— oo,— (Tg], then a'^zi9 {z)/'d{z) = a1 for all z G (0,oo). If a > —a1, then \\m.t^oo Zt = oo 
in distribution. If a < —cr^, then 

¥' [Zt G dy) -^ c y {al + aly) % dy (42) 

(weak convergence of measures on (0,oo)J for every z G (0,oo) where c = c{a,al,a'^) G (0, oo) 
is a normalizing constant such that the right-hand side is a probability distribution. 

The proof is deferred to Section [6l 

Theorem [7] exhibits a difference in the survival opportunities between the weakly subcritical 
and stronlgy subcritical regimes. In the strongly subcritical regime, ?? is a linear function and 
the SDEs of the conditioned process {Zt, St)t>Q simplify to 

dZt = (al + \(ylZ^ dt + ZtdSt + J^tdwf^ 

\ ^ J (43) 

dSt = (a + al) dt + aedwi^'^ 

for t G [0,oo). The drift of the associated Brownian motion {St)t>o is increased by al, but 
is still negative. Thus the conditioned process survives solely due to the immigration term 
aldt. In the weakly subcritical regime, the drift of the associated Brownian motion {St)t>o is 
strictly positive. Thus the conditioned process in the weakly subcritical regime survives due to 

a supercritical environment with positive probability. The immigration term al^y^Ztdt is not 
needed for this. The effect of this immigration term is to ensure survival with full probability. 
Another observation in the weakly subcritical regime is that the environment in the conditioned 
process depends on the population size. The reason for this is that the survival probability of a 
supercritical BDRE depends on the initial mass. In addition it is intuitive that the environment 
in the conditioned process needs to be less beneficial if the population size is large. Formally this 
means that the additional drift term alz-d {z)/'d{z) is decreasing in z G (0,oo). More precisely, 
this function decreases in the weakly subcritical regime from al to —a as the population size 
increases . Moreover, Theorem [7] provides an explicit quantification of the dependence of the 
additional drift term in the conditioned process on the population size. 

In the strongly subcritical and in the intermediately subcritical regimes, conditioning on ul- 
timate survival affects each individual in the same way so that the conditioned process is again 
a BDRE except for an additional immigration term. In the case of a constant environment, 
branching processes with immigration may be represented as a branching process with an addi- 
tional immortal individual. The trajectory of the immortal individual is referred to as spine or 



backbone. This backbone construction goes back to Kallenberg (1977) for branching processes 
in discrete time and has later been established e.g. for branching processes in continuous time 
(Gorostiza and Wakolbinger 1991), for the Dawson- Watanabe superprocess (Evans 1993), for the 
infinite-variance (1 + /3)-superprocess (Etheridge and Williams 2003) or for general continuous- 
state branching processes (Lambert 2007). In all of these backbone constructions, families evolve 
independently of each other. The next theorem establishes the backbone construction for BDREs 
conditioned on ultimate survival in the strongly subcritical and in the intermediately subcriti- 
cal regimes. Here the families are correlated through the environment. Due to Proposition [21 
however, this correlation is rather explicit. 

We understand a family to be a single ancestor together with the progeny of that individual. 
The total mass hereof as a function of time is an excursion from as a single individual has mass 
in the diffusion approximation. We denote the space of continuous excursions from as 

[/ := {X G C ((-oo,oo), [0,oo)) : Tq{x) G (0,oo], Xi = Vt G (-oo,0] U [To(x),oo)} (44) 

where Tq{x) := inf{i > 0: Xt = 0} G [0,oo] is the first hitting time of (inf := oo). Let {Ft)t>Q 
be the strong solution of the SDE 

dFt = ^/FtdWt (45) 

for t G [0, oo). The law of families of the process {Ft)t>o is the excursion measure Q which is a 
cr-finite measure on the excursion space U and which is uniquely determined by 



/ 9 (X) Q{dx) = lim- ^' \g m)t>o) ] (46) 



for all bounded, continuous functions g: C ([0, oo), [0, oo)) — >• IR depending on a finite time 
interval and such that there is an e > such that g{x) = for all x G C ([0,oo), [0,oo)) with 
suPi>o Xt < e. 

Theorem 8. Assume cjf,, Ge G (0, oo) and a G (— oo, — Cg]. Let {W^ )t>o be a standard Brownian 
motion. Define 

• St := [a + al)t + UeW^'"'^ for t G [0, oo), 

• a random time {f{t))t>o through f{t) := Jq e~^^a'^du for t G [0,oo), 

• a Poisson point process V on [0, oo) x \J with intensity measure dy x Q and 

• a Poisson point process V on [0, oo) x U with intensity measure dt x Q. 

Assume the ingredients {W^ )t>o, V andV to be independent. Let {Zt, St)t>o denote the BDRE 
{Zt,St)t>o started in Zq = z £ (0, oo) and conditioned to never go extinct and define a process 
{Zt)t>o through Zq := z and through 

Zt-= Yl h<zXme'^'+ Yl Xr{t)-ue^' (47) 

(y,x)&v («,x)6P 

fort G (0, oo). Then {Zt,St)t>o o-i^d {Zt,St)t>o are equal in distribution. 



10 



The proof is deferred to Section [71 

The excursions {y, x) &V are the famihes whose ancestor hved at time 0. Due to condition- 
ing on ultimate survival, there is an immortal individual. Offspring of this individual are the 
ancestors of families (s, x) G V. Conditioned on the environment, all of these families evolve inde- 
pendently of each other. Note that the environment appears in (|47p only through the time-change 
and through the reweighting of the critical excursion paths. In the critical and weakly subcritical 
regimes, the conditioned process {Zt, St)t>o is not a BDRE with immigration. A representation 
with independent families is therefore not possible. In view of Theorem [8l the SDEs ()40p can 
still be interpreted as follows: Birth events of the immortal individual are accepted only with 
probability zt!) {z)/'d{z) € (0, 1) if the current population size is z G (0, oo). Moreover the addi- 
tional drift of {St)t>o is not a^ as in the strongly subcritical regime but alz^ {z)/-&{z) G (0,o"e) 
if the current population size is z G (0, oo). 

2 Diffusion approximation 

Proof of PropositionUi We derive Proposition [T] from Corollary 2.18 of Kurtz (1978). To check 
the assumptions hereof we need more notation. Define 



a 



W ._ 



E 

fc=0 



k 



m 



iQl 



(n)) 



i)'Qr^(^) 



(48) 



(n) 

for all i G N, n G N and note that a^ , i G N, are independent and identically distributed for 
every n G N. According to assumption (fTOjl . the expectation of a^'^ converges to cr^ as n —;• oo 
for every i G N. Therefore the law of large numbers for triangular independent sequences implies 
that 

^ lnt\ 

An{t) := - ^ «!"■' ^tal as n -;> oo (49) 



n 



i=l 



almost surely for every t G [0,oo). The rescaled associated random walk converges (see 
Theorem 3) to a Brownian motion {St)t>o, that is, 



^ 
>->i 



(n) 
tn 



(St) 



t>o 



t>0 



as n — )• oo. 



(50) 



The Brownian motion {St)t>o has drift a due to assumption ([7]) and due to log(x) ~ x — 1 
for all j; in a neighbourhood of 1. Furthermore, {St)t>o has infinitesimal variance Cg due to 
assumption ([S]). Moreover, Corollary 2.18 of Kurtz (1978) requires a third moment condition 
which follows from (191) and from 



E 



n 



^ [nt\ oo 
3/2 Z^ 2-^ 



k 



i=l fc=0 '^(Q 



H-, 



^{n) 



QT\k) 



< 



sup E 

n new 



E 



k 



■k=o rn{Q\ 



(nh 



Qf\k) 



Having checked all assumptions. Proposition [H follows from Corollary 2.18 of Kurtz (1978). D 
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3 The Laplace transform and the extinction probabiUty 

Proof of Corollary It suffices to prove (fT9]) for the version (fT7|) of the BDRE due to Proposi- 
tion [2j The Laplace transform of Feller's branching diffusion {Ft)t>o satisfies 



E^ 



exp (-XFt 



exp 



i^ + i 



for t,z,X £ [0, oo), 



(51) 



(e.g., Example 26.11 of t28j). Thus we get for the Laplace transform of F^u\e^* that 



E^ 



exp(-AF,(t)e**)l(5.) 



s<t 



exp 



exp 



Ht) + jk 



/p ^ exp ( - Ss)ds + J exp{-St] 
for all t,z, X £ [0, oo) almost surely. This completes the proof. 



(52) 



D 



Proof of Corollary^ Fix t £ [0,oo) and z € [0,oo). Letting A — )• oo in formula ()19p for the 
Laplace transform and applying the dominated convergence theorem yields that 



A— >oo 



(^'0)(Zi > I {Ss)s<t) = nmJE(^'O) 1 - exp {-XZt)\ {S. 
1 — exp 



s s<t 



lim 

A— >oo 



1 — exp 



/q ^ exp {- Ss)ds + J exp{-St 



(53) 



/o ^exp{-Ss)ds^ 

almost surely. Now assume that ag > and recall that /3 = — ^. Note that 

2acj2 



as - aeW^^^ 



0<s<t 



^e 4 ^s/ o<s<i 



(54) 



.(e) 



Inserting 5^ = as + aeWg , s < t, into ([53]) . applying ([M|) and the time substitution u := ^s, 
we get that 



'■(^Zt>0\{Ss)s<t) =l-exp 



d ^ ( 

= 1 — exp - 



Jq ^ exp (—as — aeWs )ds 



V /* ^ exp -2%^s + 21^5 ds 



1 — exp 



1 — exp 



jf*/'^4,exp{2(3u + 2w!f^)du 



zat 



(55) 






24^^) 
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almost surely. Taking expectations, we arrive at 



[Zt > 0) = E 



2A 






f{za)pt„2/i^f3{a)da. 



The last step is equation (2.5) in Matsumoto and Yor (2003) which requires /3 > — 1. 



(56) 



D 



4 Asymptotic behavior of the survival probabihty 

Throughout this section, let {Zt, St)t>o be the strong solution of ([T]) and / be defined as in ([2 



4.1 General results 

Recall Al" , t G [0, oo), 7 G R, from ([20|l . First we provide sufficient conditions under which 
E[z/Aj ] and E[(1 — exp(— z/^J ))], z > 0, have the same asymptotics as t — )• 00 for 7 € R. 
We will use this for the intermediately and strongly subcritical regimes. 

Lemma 9. Let {Yt)t>i be a family of non-negative, real-valued random variables. Assume that 
there exist a function c: [l,oo) — )■ [0,oo) and a constant a G [0,oo) such that limj_>.oo ctE[yj] = a 
and limsupj_^o^ ctE[y/] = 0. Then we have that 



lim qE 

t— >oo 



1 - exp ( - Ay; 



Xa 



(57) 



for every A G [0, 00). 

Proof. Let A G [0, 00) be fixed. The upper bound follows from 1 — e~^^ < Ax, 2; > and from 
limi_>.oo qE [yj] = a. The lower bound results from 1 — e~ > Xx 
Applying this, we get that 



\2 2 

^^ for every x > 0. 



qE 



1 - exp(-Ayt) 



>qE 



Ay-^ 



ActE[yi] - ct ^ 



for every t > 0. The assumptions of the lemma then yield that 



lim inf qE 



exp ( - Ayt) > hm inf ( XctE[Yt] 



A2QE[y,2 



Xa, 



which is the lower bound in ([5 



(58) 

(59) 
D 



The Cameron-Martin-Girsanov theorem (e.g.. Theorem IV. 38. 5 in [35]) for a standard Brow- 
nian motion {Ws)s>o asserts that 



E 



h{{Ws + -/s) 



0<s<t 



E 



eMlWt - jh/2)h{{Ws)o<s<t) 



(60) 



for every 7 G R, every measurable function h: C ([0,t],lR,) — )■ [0, 00) and every t G [0, 00). The 
next lemma will allow us to deduce the intermediately and the strongly subcritical case from the 
critical and supercritical case, respectively, by changing the drift through ([60|) . 



13 



Lemma 10. Let 7 G R and let (Aj )i>o he defined as in /120\). Then we have that 



E 



2A 



(7) 



E 



(2^1^ V 



= e-(2^-2)*E 
< e^(27-2)tjg 



1 



2A 



(-(7-2)) 



1 



2A 



(7-2) 
t/2 



•E 



1 



(-(7-2)) 



2A 



(61) 
(62) 



t/2 



/or every t G (0, 00). 

Proof. Let t G (0,oo) and 7 G R be fixed. Applying ([60]) . we obtain that 



E 



1 



(7) 



■2A^ 

e~2*E 



E 



2/o*exp(2(7s + VFi'^))(is 
1 



E 



exp(2w/''^ - 22t/2) 



2/Jexp(2((7-2)s + Tyr))^^ 



L2 /o* exp (2((7 - 2)s + H^i"^) - T^i')))ds- 
Next we use the substitution u := t — s and get that 



e"2*E 



L2/o*exp(2((7-2)s + H^/!i))d. 



E 



1 



2A 



(7) 



e-2*E 



2 /J exp (2((7 -2){t-u) + Wr))du 



g-(2+2(7-2))t]g 

e-(27-2)tjE 



1 



2 /({exp (2(-(7 - 2)n + wif'^))du 
1 



(63) 



2 /J exp (2(-(7 - 2)u + Wr))du 



e-(27-2)t]E 



1 



2A 



(-(7-2)) 



which is the first claim of the lemma. 

For the second moment of l/(2Aj'^ ), we use analogous arguments to obtain that 



E 



1 



(2A^¥ 



E 



1 



E 



exp(2w/"^ - 22t/2) 



2 £ exp (2(7s + W}'>))dsy^ L (2 J^ exp (2((7 - 2)s + Wr))ds) 



rie). 



A rough estimate allows us to split the integral into two independent parts: 



E 



(2A1^))2 



< e"^*E 



e"2*E 



e"^*E 



exp{2Wl 



(e)^ 



{2f^exp{2{{j-2)s + W^''^))dsy 
1 



exp(2W/'^) 



.(e) 



2 £'' exp (2((7 - 2)s + W^D)^^ 2 /* exp (2((7 - 2)s + Wr))ds 



(e)N 



t/2' 



(64) 



2A 



(7-2) 

t/2 



•E 



2 /* exp (2((7 - 2)s + W^'> - Wt^'>))ds 
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Substituting u := t — s yields that 



E 



1 



< e"^*E 



^ g-2t-2{7-2)tjg 



2A 



(7-2) 
t/2 



•E 



exp(-2(7 - 2)t) 



2j'exp{2{{j-2){s-t) + W}l[))ds 



t/ 



2A 
1 



(7-2) 
t/2 



•E 



1 



= e^(27-2)tjg ^_^ . ^ 
[94(7-2) J 

wliich is the second claim of the lemma, 



2 /J/" exp (2(-(7 - 2)u + TyD)ds 
1 



2A 



(-(7-2)) 



(65) 



D 



4.2 The supercritical regime 

In this subsection, we will prove the result of Theorem [5] in the case of S having positive drift 
a > 0. Here we will prove for a > that, starting from z > 0, the probability of survival is 
strictly positive in the limit as t — t- 00 and for every z € [0, 00) is given by 



lim P'^(Zi > 0) = E 



t— >oo 



1 — exp 



cr. 



o. 



A — T 



1-1 + 



0-; 



ot 



2a 



(66) 



where Gy is gamma-distributed with shape parameter v G (0, 00) and scale parameter 1, that is, 

1 



P(G^ G dx) 



V{v\ 



x" e "" dx, X € (0,00). 



(67) 



The last equality in (j66p follows from the explicit formula for the Laplace transform of the 
gamma-distribution (e.g. [SZ]), that is, for every A > and v G (0, 00), E[exp(— AG^)] = tjtaT^- 
We use the following result of Dufresne (1990) (see also 



(68) 



Lemma 11. Let {Bt)t>o be a standard Brownian motion. For all a ^ and b > 0, 

exp{aBs — bs) dsj =—G2b 



where Gy is gamma- distributed with shape parameter u := ^ and scale parameter 1, i.e., has 
distribution (l671). 



Proof of the supercritical case of Theorem,\^ Let z > be fixed. As 1 — e~^ < 1, x > 0, we may 
apply the dominated convergence theorem. Using (j25|) in Corollary S] and continuity of /, we get 
that 



lim ¥\Zt > 0) = lim E 



1 — exp 



2z 



cr. 



2 i-t 



/o exp ( 



as — UeWs )ds 



E 



(T. 



1 — exp ( 7^ zG 2a 



a 



where G^ is gamma-distributed with shape parameter v := -% and scale parameter 1 



(69) 



D 
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4.3 The critical regime 

Next we study the case a = 0. Here, we will prove ([29]) . that is, if a = 0, then 



lim Vt-F^'iZt > 0) 



i— )-oo 



V2 



cr: 



log 1 + ^-Z 



crt 



(70) 



.(0) 



for every z € [0, oo). For simplicity, denote At := A^ for t > 0. The next lemma yields the 
convergence. 

Lemma 12. Let g: [0, oo) — t- [0,oo) be a Borel measurable function. Assume for some constant 
c,p £ (0,oo) that g{x) < cx^ for every x >0. Then we get that 



lim Vi-E 



t— >oo 



'^^^ 



oo ^ g-a 

g{a) • , da < oo. 



27r a 



(71) 



Proof. Let 5 be a function fulfilling the conditions of the lemma. Instead of (|27|) . we will use a 
simpler expression for the density of At for t G (0,oo). According to Theorem 4.1 of Dufresne 
(2001), the density function of l/(2^t) is given by 



i^. ^ -^ 



= — F== — 1= / exp( -a(cosh(y))^- — ) cosh(y)cos(— -) dyda 
--■.pt{a)da 



(72) 



on (0,00) for every t E (0,oo). As for every x > 0, |e^ < cosh(x) < e^, |cos(x)| < 1 and 
g{x) < cx^, we obtain that 

e St g{a) 



sup 
Jo t>i 



Vty/a 



exp y - a(cosh(y))^ - — j cosh(y) cos (^— j 



dyda 



< I I ce » a'' 2 exp 
/o Jo 



ae 



•iy 



exp(y) dy da. 



Using the substitution y := ax, we see that 



1 



1 



yp-,x-2-Pe~y- dy = ^— r {p + 1) 



/o JO -^ xa+P 

for every x > 0. Now applying Fubini's theorem and ()74p yields that 



(73) 



(74) 



a^ 2 exp 



JO 

00 



exp{y) dy da = a^ 2 exp f — a— — j daexp(2/) dy 



r(p + 







.2y 



eydy = T{p+ i) 2^P 



1 \ o2p+l 



-2pj/ 



dy 



r(p + ^)2^-+^^<oo. 



(75) 
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Thus the integrals in ()73p are finite. Applying the dominated convergence theorem, we get that 



lim \/t • E 

i— >oo 



1 1 /"°° 

di^r) =limVt/ g{a)pt{a)da 



t">00 







oo />oo 

g{a) / lim ^ — 
Jo *^°° V7r2 Vo 



-v/2e8t 1 ^ ^ w N\2 y^\ 



exp I — a I cos 



My))^ 



2^ j cosh(y) cos (^— ) dy da 



ny 



/•oo /"OO /rj -I 

/ ^(a) / ^ exp I — o(cosh(y)) ) cosh(y) dy da. 

Jo Jo ^ V« ^ ^ 



(76) 



To complete the lemma, we will simplify the inner integral in the above equation. Noting that 
(cosh(y))^ = 1 + (sinh(y))^, y € M, and using the substitution x := y/2a sinh.{y) yields that 



V2 1 



exp I — a I cos 



71" Va 

°° V2 1 



_ exp [ — a 

vr V« ^ 2 / V2a 



cosh{y) dy= / ^ exp ( - a(l + (sinh(y))^)) cosh(y)dy 

Jo "^ v^ 



dx = / e '2 dx 

vra Jo 



vra 2 ^/2^T a 



for all a G (0, oo). 



D 



Proof of the critical case of Theorem\^ Fix z > 0. As f(zx) < -Z^x for every x > 0, the 
conditions of Lemma [12] are met with g = f ■ Applying Corollary U and Lemma [T2l we get that 



lim \/i • P^fZt > 0) = — lim \ j-^ ■ E 

t— >oo (Tp t—^oo V 4 



2^ 



tai/i 



f{za)^—da. (77) 
<7e Jo V 27r a 



Applying Fubini's theorem, we get for every c G [0, oo) that 



(1-e-^^) dx 

^ 



e-^^e"^ dy dx = / / e-^'^+i^^dxdy 
./o Jo Jo 

— — dy = log(l + c). 

y + 1 

Using this result with c= ^z and recalling the definition ()24p of /, we get that 



t— >oo 



0-, 



e JO 



27r a 



\/vrcre 



which proves Theorem [5] in the case a = 0. 



(78) 



lim Vt-P"(Zi>0) = — / (l-exp(-§za))^=--da = ^^log(l + ^-z), (79) 



D 



4.4 The weakly subcritical regime 

Next, we turn to the case —a1 < a < 0. We prove ([30]) . that is, if —a1 < q < 0, then 



limt3/2el^*.p-(^^>0) = 4 / 
t^oo ' (yf Jo 



f(za)(j)ji{a)da 



(80) 



for every z G [0, oo) where (/)/3 is defined in ()33p and where /? = —2a/a1. The following lemma 
yields the convergence. 
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Lemma 13. Let 7 > and let g : [0, 00) — t- [0, 00) be a Borel measurable function. Assume that 
for every x >Q, g{x) < cx^ for some b > 7/2 and c > 0. Then we get that 



lim t3/2e7^t/2]E 



t— >oo 



^(rTw) 



2A 



g{a)(l).y{a)da. 



(81) 



Proof. Let 7 > be fixed and let 5 be a function fulfilling the conditions of the lemma with 
constants b > 7/2 and c > 0. Recall the density of — j^, i > 0, from ()27p . With the substitution 

u = as, we get for t G (0, 00) that 



^3/2g7^t/2jg 



5(77^) 



2A 



(7)' 



^e^^/7 + 2 
\/27r2 



^g^V2t ^^ + 2 



00 roo roo 



JO Jo 



<7(a)e-'^a-(^+i)/2e-«'/2*f-V 



n\(7-i)/2 
a) 



\/27r2 



00 /'OO /'OO 



JO JO 



sinh(^) cosh(^) sin(7r,^/t) dw 
(n/a +(cosh(e))2)(7+2)/2 ^T " 

,-a^-7/2g-«V2t^(7-l)/2g-« 

sinh(,^) cosh(,^) sin^-K^/t) 



{u + a(cosh(e))2)(7+2)/2 



d^duda. 



(82) 



In order to apply the dominated convergence theorem, we will show that the integrand on the 
right-hand side is dominated for t G [l,oo) by an integrable function. As for x > 0, | sin(x)| < x, 
sinh(x) < cosh(x) < e^, cosh(a;) > e^/2 and by assumption, g{x) < ex , there is a constant 
d = d^ > such that 



00 POO fOO 



sup 
Jo Jo t>l 



te-V2t 



V2' 



vr^ 



r(l±l')5(a)e-«a-^/2.-«V2t„(7-i)/2,-^. 



<d 



<d 



00 fOO fOO 



e -''a^-^/2n(^-i)/2e-" 



sinh(^) cosh(,^) sin(7r^/i) 
(n + a(cosh(e))2)(7+2)/2 



JO JO 
00 />oo />oo 



(4u + oe2«)(7+2)/2 



d^ du da 
d^ du da 



e "a^'^/^n(^-^)/^e'" ?,. , ^:,„ d^duda 



-57£ 



1 JO JO 

1 POO fOO 



+ 



JO JO 



^fe-7/2^(7-l)/2g-« 



a(7+2)/2" 



(4Me-2« + a)(7+2)/2 



d^ du da 



(83) 



The first summand on the right-hand side is finite as ^^g" > — 1 and 7 > 0. Recall b — 7/2 > 0. 
Choose < e < min(6 — 7/2, 7/2). Using a^ < ay for every a G [0, 1] and < y < x, and applying 
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Fubini's theorem, we estimate the second summand on the right-hand side of (|83p as foUows: 
"1 roo roo e~^'^£ 



^6-7/2^(7-l)/2g-« 



^0 ^0 

rl /"oo /"oo 

-/ / / 

Jo Jo Jo 



^.^(7-l)/2g-n 



(4ue-2« + a)(7+2)/2 



(4ue-2C + a)(7+2)/ 



d^ dn da 
fi,^ du da 



2Ci„^S,(7-l)/2g-« 



< / / / (4^6-^^ + a)'u 
Jo ^0 ^0 

JO JO Jo 



-C7^ 



(4Me-2« + a)(7+2)/2 
1 



da du d^ 



(4ne-2€ + a)(7+2)/2- 



da du d^ 



u 



< 



7/2 - e Jo Jo 

^-7/2+e /-oo /-oo 

7/2 - e X io 



(7-l)/2g-«g-C7^ 



1 



1 



(4ne-2€)7/2-. (4ue-2« + 1)^/2- 



dud^ 



(84) 



4— >oo 



Thus, applying dominated convergence in ([82]) and using tsin(7r^/i) — > vr^ for every ^ G [0,oo), 
we get that 



hm t3/2e7^t/2E 



1 



^2A 



(7)- 



r 



7 + 2 



00 /"OO /"OO 



lini l£^-g(a)e-"a-^/2e-« /2*^(7-i)/2g-« 



JO Jo *-^°° \/27r2 

sinh(^) cosh(^) sin(7r^/t) 



1 /^ _(_ 2^ '■°° '■°° '■°° 

r' ' 



v2vr ^ 2 / JO JO JO 
g{a)(j)^{a) da, 



(n + a(cosh(0)2)(T+2)/2 



5(a)e-'^a-^/2^(^-^)/2e 



d,^ du da 

,^n sinh(^) cosh(^)g 
' (n + a(cosh(e))2)(7+2)/2 



dS, du da 







which is the claim of the lemma. 



(85) 
D 



Proof of Theorem in the weakly subcritical case. Assume — o"2 < a < 0. Let z > be fixed and 
(3 = —2a/a^. By Corollary HI we obtain that 



lim t^/^e 



t— >oo 



^ • F^Zt > = ^ lim ^^ ^^^^ ' ' 
ai t-S'oo \ 4 



2 -^2^ 



/ 



2A 



{/?) 

toi/4 



(86) 



As for every x > 0, f{zx) < -^x and < /3 < 2, the assumptions of Lemma [13] with g = f, 
7 = /3 and 6 = 1 > ^ are fulfilled. Applying Lemma [13] to ()86p proves Theorem [5] in the weakly 



subcritical case. 



D 



The following lemma for the weakly subcritical case will be needed later. Recall /3 = —2a/a1 
and 'd from (1391). 
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Lemma 14. Assume at, (T(, € {0,oo) and a ^ {—a'^,0). T/ien ■!? G C°° ([0,oo), R) and 



lim 



^{z) 



'^^zhogiz) /5 (jS sin (vrf ) V^^i 



cr. 



^\ 2 



(87) 



lim 



z^'{z) 






Z2log{z) cr| sm(7r^ 
Proof. In order to prove "i? G C°° ([0, oo), R), note that 



2\ 2 



5"Z 



a) da 



oo 2n 

6 T? 

CTfe 



^2 \ ^2n |.oo 

" za](t)fi{a)da< -^ j a"^ (t)[s{a)da (89) 



(T 



CT; 







for all z G [0, oo) and all n G N. The right-hand side of (j89p is finite according to Lemma [T3] 
(note 1 > 13/2). The dominated convergence theorem thus implies that i? G C°° ([0,oo), R). 
Moreover, we get that 

C /-oo 

(90) 



^W(z) = ^/ /(")(za)o"(/'^(a)da 



for all z G [0, oo), ?i G Nq. For proving (f87|l . we substitute b := za and get that 
i}{z) 1 8 



Z2 lOg(z) 2:2 log(z) "e JO 



f{za)(l)fs{a)da 



f^e JO 



/(b)- 



Z2^ 



log(z) V^ 



for all z G (0,oo). Next we rewrite the definition ()33p of (j)^ and see that 



db 



(91) 



z2+J^log(z) 



V2: 



' r(^±^]e-h-^/^ 



vr 



00 fOO 

^(/3-l)/2g-« 



sinh(^) cosh(0^ , . , 

2 /3+2 "s du 

zlog(z) (« + ^(cosh(^)) )^~ 



(92) 



for all b G (0, 00) and all z G (0, 00). Using the substitution x := cosh{S,) \/b/{zu) and noting 
that dx/d^ = sinh(^) y^b/{zu), we obtain that 



sinh(^) cosh(^)^ 



-d^ 



zlog(z) (n + j(cosh(^)) )^~ bu"^ Jo log(2;) (l + 



cosh(OA/^-e-sinh(OA/^ 



|,(cosh(C))')^ 



de 



°° xarcosh (xy'^^j 



^ / n~ /9+2 

b^2 ^V^log(z)(l + x2) — 



dx 



(93) 



for all b G (0, 00), u G (0, 00) and all z G (0,oo). Note that arcosh(y) = log(y + y^y^ — 1) for 
all y G [1,00) and therefore arcosh(y) < log(2y) for all y G [l,oo). Consequently, we have that 
arcosh(xyf )/log(^) < (log(xyf) V 1) + log(2) + ig^ < 3(log(xyf) V 1) for ah z G [3, 00) 
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and all x G {y^b/{zu), cxo). Let us prove that we may apply the dominated convergence theorem. 
From dUD, (PD and ([M]) we see that with c := jj-^'^m^) e (0,oo), 

X arcosh (x \/^) 



Z2 log(z) 



f{b)e~-b~^' 







poo 

Jo 



u 2e 



n 2e " 



'^log(z)(l + 2;2)^ 



dx du db 



13+2 
(1 + X2) — 



3 log x./]i VI Uxdudfe 



for ah z E [3, oo). The right-hand side is finite as /3/2 G (0, 1) and as f{b) < {bal/al) A 1 for 
b G [0,oo). Thus we may apply the dominated convergence theorem and get that 



lim 

^^~ Z2 log(z) 

j-OO 

= c f{b)r 

Jo 

/•oo 

= c f{b)b- 
Jo 



u 2e 



lim e 2 1 



arcosh (xy^) 



X 







<x 



db- 



u '^e ^ du- I - 
Jo 2 



{l + x 



8+2 



log(z) 
■ dx 



(l + x 



i^ 2 



dxdudb (94) 



For the last step we used that arcosh(z)/log(z) — )■ 1 as 2; — )• 00. We will simplify the integrals in 
the above equation. Inserting c = ■^-h-^{2~) yields that 



^(z) 



hm „ _ 

^^°° Z2 log(z) f^e V2: 



U(¥)f/Wr«-'<i..r(i).i5 



(95) 



Using integration by parts, we obtain that 

/■oo . 

f{b)b~^/^-Ub= / 6-^/2-1 fl-exp 
Jo ^ 



o-; 



CT: 



1 



-/3/2 



6-^/2 fl 



exp 



crz 



crt 



db 



2a| 



/ 61-/3/2-1 exp (-^6 




(i6 



/ —— 6-/^/2 ^exD 



cr: 



at 



6 M6 



(96) 



2ai 



/3a: 



,-2 \ ^ 2 



r(i-f). 



where we used the integral in the gamma-distribution with shape parameter 1 — f ^^'^ scale 
parameter ^. Inserting ([96]) into ([95]) and using r(^) = ^/^T, T{1 + f ) = f r(f ) and Euler's 
reflection principle r(x)r(l — x) = 7r/sin(7rx) for x G (0, 1), we get that 



lim 



^z) 



1 /3p//3 



z— ^-oo ^^ 



2a„^ 



0-; 



9N 1-^ 
2 \ ^ 2 



r 1 



/3 



z2\og{z) 0"f\/27r2 V2/ /3(t2 Va2 
This proves (j87p . Paralleling the above arguments yields 



2V27r 



a: 



2 7 V2;2/3 c73/3sin(7rf) V^6 



2\ 2 



D 
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The following lemma implies that the additional drift term in ()40p is strictly decreasing in 
the weakly subcritical regime. 

Lemma 15. Assume (Tb,(Te G (0, oo) and a G (— <Tg,oo). Then the function (0, oo) B z ^-^ 
a'lz'd {z)/i}{z) G IR, is strictly monotonic decreasing and satisfies \\Ta.z~>Qa1z'd [z)/'d{z) = a^ and 
lim2_5.oo 0"gZi? {z)/'d{z) =max(— a,0). 

Proof. In all regimes a G R, we have that •& (0) G (0, oo) and t?(0) = and, therefore, 

In the supercritical regime a > 0, we see that 
..'(.) «i(l-0 + S-)"*) ^t(l + f-)"*" 2a. 



l-(l + f|.z) - l-(l + ||.z) - ' ^ (l + f|..)--l 

for all z G (0,oo). The derivative hereof is strictly negative 

d z-d' {z) _ {ol + alz) 2a - 2azal 1 2az ~%\}^^'^, 



+ 



171 






2q 1 , , 2a 



i „ __ / _2 \ ^7 „ / _2 \ — T — i 



/o^4a2(i + ||.y)^- dy_4a2;o^(l + f|-2/)'^+y!f(l + |l-y)'^~ ^2/ 



+ -e^^)'((l + fl-)^-l)' 



2q 1 
t2 \T7-1 



- -^ 2 <0 



for all z G (0,oo). Moreover, it is clear that lim^_!.oo zi} {z)/'d{z) = 0. 
In the critical regime a = 0, we get that 



z^\z) -i^(l°g(l + f|-)) al 



azz 



^(^) log f 1 + 42) ^b + "^^^^ log f 1 + 4^ 



6 / \ "^b 
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for all z £ (0, oo). The derivative hereof is strictly negative 

d zi^' (z) _ {a^ + a^z) al - a^zal 1 ^ a^z -1 a^ 



dz m {al + alzf log (l + ||.) ^^l + ^l^ (log (l + ||,))' ^' + '^l^ 



aWAo^il + ^z] -alalz 



jOglUg ^i-h -^z\ - o^o^. 

^ < 



(<r2 + alif (log (l + 5 



*^" 



for all z G (0, oo) where we used the inequality log(l + x) < x for all x E (0, oo). Moreover, it is 
clear that lim2_^oo zd {z)/'i}{z) = 0. 

For the rest of the proof, we assume that a G (— (Tg,0). Recall /3 = —2a/a'^. Lemma [HI 
implies that 



2 Z'& {z) 2 2^ log(^) 2 '^g'^i"(^|) '^'^"^ 2/3 /oe^ 



.^log(.) /3,3,i,(,|) V^ 



It remains to prove monotonicity of (0, oo) 9 z i— )■ 21? {z)/'&{z). The first derivative hereof is 
d zd\z) _ ^(^) (^» + ^^'V)) - z^'{z) ■ ^'{z) 

(f T?(Z) - Zl?' (Z)) ^'{Z) + ^(Z) ((1 - f ) ^'(^) + Z^"{Z)^ 

for all z € (0,oo). We will show that the right-hand side is negative for all z G (0, 00). Define 
a function (^^ : (0, 00) — )■ [0, 00) through (^/^(a) := -T-(/)^(a)a2+^ for all a € (0,oo). Using the 
substitution r] = ^/a cosh{S,) , we rewrite this function as 

a^7o Vo V27r V 2 y (u + a(cosh(e))2) 2 



2tt y 2 J Jo J^ (u + n'2)^ V«sinh(0 



cT^\/27r \ 2 J Jq J^ (u + r?2) 

' ' :r(^)e-"/°°^^^-^^^^^-"/°°Vv^arcosh(-^) '....2 ^^^^ 



cjlV2n \ 2 J J, 7o ''^^'^ V^/^y (n + r?2)^ 

for all a G (0,oo). As (0,oo) 9 a 1— )• e~" arcosh (-?=) 1^ zjj is strictly monotonic decreasing for 
every ij G (0, 00), we conclude that 0^ is strictly monotonic decreasing and that (j)g{a) < for all 
a G (0, 00). Moreover (j)i3{a) decays at least exponentially fast as a — )• 00 so that / lP(j)p{h) dh < 
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oo for all a,p (^ (0,cx3). Applying ()90p and integration by parts twice, it follows that 

o /-oo /-oo 



iz) = — / /(")(za)a"(^fl(a)(ia 



oo 



(•(")(' — ^ / un-i-l. 



r\za) / 6"-t-V/3(6)rf6 



(n+i)/ . / ro 2 



/^"+^^(za)z 



n 



<^M&) 



oo poo poo 

+ / /("+i)(za)z / V'^2-^^p{h)dhda 

JO Ja 

-g4>p{h) db da 



o- Ja n 



n — IT JO 



f^'^+^\za)za'^^^^p{a)da-z / Z^+^^M 



■g4'p{h) dbda 

a n— TT 



(100) 



for all z G (0, oo) and all n G No- Recall that 0o(6) < for all h G (0, oo) and note that 
f'^'^\za)^^p{b) > for all a,b,z £ (0,oo). Thus (fTOOjl with n = implies that 



f.«<-f^f- 



-^ 7 



OO 



f^'->{za)za 2(j)i^(^a)da = z—^ p>{za)a(j)p{a)da = z'd{z) (IQI) 

<^e iO 



for all z G (0,oo). Moreover /(2)(za)^(^l(6) > for all a,b,z e (0,oo). So ([TOO]) with n = 1 

2 

yields that 

l-f)/(z)<- / f^-^\za)za^-2 4>f^{a)da = -z^ / /(2)(za)a2</.;3(a) da = -z/'(z) 

^ Jo ^e Jo 

(102) 



for all z G (0, oo). Applying the inequalities (jlOip and ()102p to the right-hand side of ([99]) and 
using -d {z),i!}{z) > for aU z G (0,oo), we conclude that ■j^{z'& {z)/'&{z)) < for all z G (0, oo) 
which implies that (0, oo) B z >-^ ztl) {z)/'d{z) is strictly monotonic decreasing. D 



4.5 The intermediately subcritical regime 

In this subsection, we will prove ([3T]) . that is, if a = —a1 (/3 = 2), then 

V2ae 



limVte^^F''{Zt > 0) 



t— >oo 



vrcr 



2' 



(103) 



for every z G [0, oo). 



Proof of Theorem\^in the intermediately subcritical case. Let z > be fixed for the moment. 
Lemma do] asserts that ^[^py] = e-2*E[2^] and E[ — 5__] < e-2t(|]g^_^j^2 ^^^ ^ ^ (o,oo). 



L2A 



(2AJ"')^- 



-lt/2J 



By Lemma [121 limj_j.oo V^^^iz^] ~ /o°° ^~k"^~^ da < oo and, therefore, (^[2^^^ — ]) decays like 
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t ^ as t — 7- oo. Thus the assumptions of Lemma [9] are met with q = \/te^* and Yt = l/(2Al ), 
t > 1. Applying Corollary U] and Lemma El we get that 



lim vie 2 

t— >oo 



{Zt > 0) 



J 


4 


lim 


>f.^# 


E 


V 


<^i 


t— >oo 


V 4 




2 




-e^ 


/■°° 1 e 


—a 


0"e 




-.^^ 


\/2^ 


a 



/ 



2yl 



da 



crt 



which proves Theorem [5] in the case a = —a^ 



(2) 
tai/4 

2a, 1 



27r' 



(104) 
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4.6 The strongly subcritical regime 

Finally, we will prove ([32]) . that is, if a < — Cg (/3 > 2), then 



lim e-("+^)*-P^(Zt >0) = ^2- 



t— >oo 



-a — a: 



0"t 



(105) 



for every z £ [0,oo). 



Proof of Theoreml^ in the strongly subcritical case. Let /3 = —2a/a'^ and assume a < — (Jg. Let 
t > and z G [0, 00) be fixed. The main tool for the proof is Lemma El Let us check the 
conditions of that lemma. Using Lemma [TOl we obtain for the first and second moment of 
l/{2Af^) that 



E 



2A 



m 



E 



L(2^ 



(/3)^2 






1 



("(/3-2)) 



2A 



1 



2A 



(/3-2) 
t/2 



•E 



1 



(-(/3-2)) 



L2A 



(106) 
(107) 



t/2 



for every t G (0, 00). As /3 > 2, the monotone convergence theorem and Lemma [TTI yield that 



hm E 

i— >oo 



1 



2A 



(-(/3-2)) 
t/2 



E 






E 



G 



/3-2 



/3-2 <oo 



(108) 



where Gp-2 is gamma-distributed with shape parameter /? — 2 > and scale parameter 1. Using 
relation (1.1) from |33] and monotone convergence, we get that 



hm E 



1 



^2 A 



(/3-2) 



lim E 



1 



^2 A 



(-(/3-2)) 



E G/3„2 = 0. (109) 

Thus, by (jl06|) . (jl08|) and (jl09|) . the assumptions of Lemma O are met with q := 6^2/3-2)* g^j^^ 
It = l/{2Af^), i > 1. By CoroUaryU LemmaEland (HUH]), we obtain that 



lim e(2/3-2)-e*/4 . p-(Zt > 0) = lim e(2/3-2)-e*/4 • eT/ 



t— ^00 



t— >oo 



2A 



2;(T" 



hm E 



at t^oo 



2A_2 



(-(/3-2)) 
t<72/4 



(/3) 
ta2/4 

4(/3-2) = .2--"-^' 



(110) 



crt 



err 
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Note that (2/3 - 2)o-2/4 = -(a + ^). Thus, (fTTOl) is the claim of Theorem [5] for the strongly 



subcritical case. 



D 



5 Proof of Lemma [6] 

Proof of Lemma\^ Fix t G [0, oo) and x £ I for the moment. As 'r]{x) > 0, there exists a 
To e [0, oo) such that P^(Xt ^ A) > for all T G [Tq, oo). Let (^: C ([0, t], /) ^ IR, be a bounded 
and Borel measurable function. The Markov property of {Xs)s>o implies that 



E^[cPi{X,)s<t)lXr^,^A] 



lXr+,^A 



(X. 



= E^'[0((X,),<i)P^*(XT0^) 

for all T € [0,oo). Consequently, we get for the conditional expectation that 

E-[</)((X,),<t)P^'(Xr0A)] 



s s<t 



(111) 



E^ 



^{{X. 



s)s<t) 



Xr+t A 



P- {Xr+t A) 
W [,/. ((X,),<i) (?(r)e^^P^* (Xt A)] g(T + i)e^(^+*) 



(112) 



q{T + t)e^(^+*)P^ (Xr+t ^) g(r)e^^ 

for all T G [T^oo). Due to suprg[i o^) q{T)e^^W^' {Xt ^ A) < c{l+ \\Xt\\P) for some constant 
c G [0, oo) and due to E[||Xi||P] < oo, we may apply the dominated convergence theorem and 
obtain that 



E^ 



limE^ 

T-5.00 



^{{X^ 



s s<t 



Xr+t A 



iiXs)s<t) lim g(r)e^^P^* {Xt A) 



lim g(r + t)e>^(T+t)-px ^Xr+t ^) 

T— s>oo 

K^[<P{{Xs)s<t)v{Xt)] 



■ hm , ' e 

T^oo g(r) 



-At 



ri{x) 

li (j) = 1, then the left-hand side is equal to 1 and, therefore, the right-hand side is equal to 
1. Consequently the right-hand side defines a probability distribution on C ([0, i],lR, ) for every 
t £ [0, oo). These probability distributions are consistent for t G [0, oo). So Kolmogorov's 
extension theorem (e.g. [28]) implies existence of a stochastic process {Xt)t>o having continuous 
sample paths and being uniquely determined by 

E^[<Pi{Xs)s<t)viXt)] 



lim E^ 

T-s>oo 



4> {{Xs)s<t) Xr+t ^A =E^[^ {{Xs)s<t)] 



-xt 



r]{x) 



(113) 



for all t G [0,oo). This proves ()36p . Moreover if t G [0, oo) and if g. I — ^ [0,oo) is a Borel 
measurable function, then (j38|) follows from (jllSp with (j){{Xs)s<t) '■= min ((7(Xt),n), n G N, 
and from the monotone convergence theorem as n — )• oo. 

Next we identify the linear operator of the martingale problem solved by {Xt)t>o- Similar 
arguments as above imply that 

rj{x) = = lim q{T + t)e^(^+*)p^ (Xr+t A) 

(114) 
lim g(r)e^^P^* {Xt A) 



lim 



^iT + t)^xt^. 



T->oo q{T) 



T-^oo 



e^'E^ [7]{Xt)] 
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E^ 



for all X G / \ / and for all t € [0,cx3). Moreover, the Markov property of {Xt)t>o, the relation 
E^[7/(Xt)] = e~^^rj{x) for x G / and t G [0, oo) and Proposition 4.1.7 of [l7| imply that 

dr]{Xt) = -\r]{Xt) dt + dMt (115) 

for t G [0,oo) where (Mj)t>o is a suitable martingale. Let (7: / — t- IR, be bounded and twice 
continuously differentiable. Ito's lemma shows that 

dg{Xt) = (Gg) (Xt) dt + (Vga) (Xt) dWt (116) 

for all t G [0,00) where Qg := Vg fi + ^ (Vcru^V*) g. Thus Ito's lemma and symmetry of aa^ 
result in 

de^'r]{Xt)g{Xt) 

= Xe^'rj{Xt)g{Xt) dt + e^* dr]{Xt)g{Xt) 

= \e^''n{Xt)g{Xt) dt + e^'r^iXt) dg{Xt) + e^'g{Xt) dr,{Xt) + e^* {VT^aa'V'g) {Xt) dt 

= e^'iiiXt) {Gg) {Xt) dt + e^* {r]Vga) {Xt) dWt + e^'g{Xt) dMt + e^* {Vr]aa'V'g) (Xt) dt 

for all t G [0,00). Taking expectations, we infer that 

e^'7]iXt)g{Xt)] - v{x)g{x) = / E" [e^"r?(X„) {Qg) {X^)] du (117) 

for all f G [0,00) and all x G / where Qg := Qg + -V7]aa'''V^g. This implies for the Markov 
process {Xt)t>o that 

E- [5(X0] - g{x) = / E- [(^5) (X)] rfn (118) 

Jo 

for all t G [0,00) and all x £ I. Now Proposition 4.1.7 in [T7] implies that {Xt)t>o is a solution 

of the martingale problem for Q. Finally, Theorem V.20.1 of |36] shows that {Xt)t>Q is a weak 

solution of the SDE (j37p . This completes the proof. D 

6 The BDRE conditioned to never go extinct 

Proof of Theorem^ Fix a G R and o'bi'^e G (0, 00). We will prove Theorem [7] by applying 
Lemma [6] to the process {Xt)t>o = {Zt,St)t>o which has state space / := [0,oo) x R. Define 
/i: /-;> 1R2 and ct: I ^ 1R2x2 i^y 

/.(.,.) = ( («+]^')^ \ and a{z,s) = ( ^' ^^A (119) 

for all {z,s) G [0,oo) x R. We set ^ = {0} x M. Note that {{Zt,St) ^ A} = {Zt > 0} for all 
t G [0,00). Moreover, define r]: I ^- [0,cx3) through r]{z,s) := i9{z) for all {z,s) G [0, 00) x R. We 
will check the assumptions of Lemma [6] for the different regimes separately. In all cases we have 
that 

E"[|Zt|2] + E^[|5tp] < ^^E^ie^^*] + 2s^ + 2E°[S^] = z^e^'''+'''' + 2s^ + la^t + 2aH^ < 00 (120) 

for ah {z,s) G [0,oo) x R and ah t G [0,oo). 
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The supercritical regime Let g = lGQ, A = and p = 0. Theorem [5] implies that 



lim g(t)e^* p(^'")(Zi > 0) = 7]{z, s) = ^z) = 1 - fl + 



t—>-oo 



cr; 



crt 



■ z 



2a 



fl21) 



for all {z, s) G [0, oo) x R. The function rj is twice continuously differentiable and satisfies 
r]{z, s) > if and only if (z, s) G (0, oo) x R. Moreover, it is clear that g(i)e'^*P(^'*) {Zt > 0) < 1 
for all {z, s) £ I and for all t G [0, oo). 

The critical regime Define q £ Q through q{t) = \/t for t G [0, oo) and let A = and p = 1. 
Theorem [5] implies that 



lim q{t)e^^F^'''\Zt > 0) = r]{z,s) = ^(z) = -£- log (l + ^-z 



t—^oo 



\/vro-e 



(Tt 



(122) 



for all {z, s) G [0, oo) x IR and thus r] is twice continuously differentiable. Note that r]{z, s) > 
if and only if {z,s) G (0,oo) x R. Moreover, Corollary [3] implies that 



1+ U,s 



sup q 



'*JII te[i,oo) 



(^t)e^tpiz,s) ^^^ > Q) 



1 



l + \\iz,s 
z 



sup 

'•^Vll i6[l,oo) 
2 



VtE 



f 



2A 



tai/4 



< sup Vt -| E 



l + ^te[l,oo) O-fe 1-2^4^2/4 



(123) 



for all {z,s) G [0,oo) x R. The right-hand side is finite according to Lemma[T2]and is uniformly 
bounded in (z, s) G [0, oo) x R. 

The weakly subcritical regime Define g G Q through q{t) = y/t for t G [0, oo) and let 



and p = 1. Theorem [5] implies that 



lim q{t)e^^ F'-'''\Zt > 0) = ??(z, s) = ^{z) 



3 / f{za)(l)p{a)da 

'^e Jo 



(124) 



for all {z, s) G [0, oo) x R. The function rj is twice continuously differentiable according to 
Lemma [T^ Note that rj{z, s) > if and only if (z, s) G (0, oo) x R. Moreover, Corollary U] implies 
that 



sup q 



1 + 11(2^,5)11 te[l,oo) 



(t)e^*p(^'") {Zt > 0) 



l + ll(2>'5)ll te[i,oo) 



sup Vi e^ E 



/ 



2A._2 



(/3) 



< 



(JZ 



1 + za'^ a? 



sup 

e <^b tG[l,oo) 



t^l ('i^)-'4-i 



e "<= 



E 



1 






for all {z,s) G [0,oo) x R. The right-hand side is finite according to Lemma [T3l and is uniformly 
bounded in {z, s) G [0, oo) x R. 
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The intermediately subcritical regime Define q ^ Q through q{t) = y/i for t G [0, oo) and 

2 

let A = ^ and p = 1. Theorem [S] implies that 



lim g(i)e^*p(^'')(Zt > 0) = r]{z,s) = t?(z) = z 



V2, 



(Te 



t— >oo 



V^af 



(125) 



for all {z, s) G [0, oo) x R. The function rj is twice continuously differentiable and satisfies 
ri{z,s) > if and only if {z,s) £ (0,oo) x R. Corollary H] implies that 



1+ U,s 



■ sup g(t)e^*p(^'^) {Zt > 0) 
te[i,oo) 



1 + IK2,S)|| t6[i,oo) 



sup Vi e a" *E 



/ 



2A 



(2) 



< 



" '— sup x/^e^^E 

b \ ^e tG[l,oo) 



1 + z 0-? V 0-2 ^c^^'!^^ V 4 



ta2 



e „2 



trTp 



2A 



(2) 
t<72/4 



for all {z, s) G [0, oo) x R. The right-hand side is finite according to Lemma [TD] with 7 = 2 and 
according to Lemma [T2l and is uniformly bounded in {z,s) G [0,oo) x H. 



The strongly subcritical regime Let q^ = 1 G Q, A 

Theorem [5] implies that 



a + 



(2/3-2)^ andp = 1. 



lim g(t)e^*p(^'')(Zt > 0) = r]{z,s) = t?(z) = z • 2- 



-a — ar 



t— >oo 



(126) 



for all (z, s) G [0, 00) x R. The function rj is twice continuously differentiable and satisfies 
r]{z,s) > if and only if {z,s) G (0, 00) x R. Corollary U] implies that 



1^11! ,11 sup g(t)e^¥(^'^) (Zt > 0) = J ™, 

l + ll(2>s)ll te[i,oo) 1 + 11(^,^)11 te[i,oo) 



sup e^*E 



/ 



2A 






<^4 sup e(2/^-2)^E 

1 + ^ O-j t6[l,oo) 



2a:_2 



toi/4 



2; (t: 



1 + ^ <^6 tG 



sup E 



6 tG[l,oo) 



2A_2 



(-(/3-2)) 
ta2/4 



(127) 



for all {z, s) G [0, 00) x R. The last step follows from Lemma [TOl The right-hand side of (|127p is 
finite due to — (/3 — 2) < and due to Lemma [TT] and is uniformly bounded in (z, s) G [0, 00) x R. 



Application of Lemma [B] After having checked all assumptions, we apply Lemma EJ The 
additional drift term is 



'cra*V*r/) {z,s) 



ri{z,s 



1 

W) 



'alz 
aeZ ae 







alz + alz"^ 
alz 



(128) 



^{z) 
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for {z,s) £ I = (0,00) X R. Inserting this into (f37|) . we get for {Zt,St)t>o that 

dZt = ^^ (alZt + alzf) dt + Qa^Zj + aZ^ dt + .pl^tdwf^ + Zta^dwf^ 

dSt = ^-^alZt dt + adt + aedW^''\ 

Therefore {Zt,St)t>o solves the SDEs ()40p . Moreover, Lemma [6] impHes that the conditioned 
process satisfies (j4ip . In addition Lemma [T5l establishes the properties of the function (0, 00) B 
z^alzd\z)/'d{z). 

It remains to establish the limit of Zj as t — )• 00. Note that {Zt)t>o is a one-dimensional 

diffusion with drift term /i(z) := -/u^io'lz + 0"gZ^) + ^a^z + az, z G (0, 00) and diffusion term 
a'^{z) := a^z + a'^z'^, z G (0, 00). Define a scale function R: [0, 00] — )• [—00,00] through 

R{z) := ^' exp (- J' ^ dz) dy (129) 

for all z € [0,oo]. Standard results (e.g. [26]) show that Z^ — )• 00 in distribution as t — )■ 00 
if -R(O) = —00 and R{oo) < 00. Let a G R. We rewrite the integral in the exponent on the 
right-hand side of (jl29|) as 

y2^{a^hZ + a^,z^)+a^,z + 2az H ^' (z) , fyal + 2a 



dz = 2 — f-^ dz + / -^ ^ dz 



r 

Ji (^Iz + crlz"^ Ji ^{z) Ji crl+a^z (^gg) 

= log {(mf) + ^^ log (^.^ + a'eV) - log ((^(1))^) - ^^ log (a,2 + a!) 
for all y G (0, 00). By ()130p . there exists a constant c G (0, 00) such that 

^(^) = ^ /' 7^7^ (^' + ^e2/)"^ ^y (131) 

Ji (??(y)) 

for all z G [0, 00]. As ■d{z)/z — )• c as z — t- for a constant c = c(a,cre,crft) > 0, we have that 
lim5;_j.o i?(z) ~ /^ -^ dy = —00. Next we show that R{oo) < 00 whenever a > —al- If a > 0, 

then •d{z) = 1— (l + %-z) "^^ and thus R{oo) < 00. In the case a = (and thus /? = 0), we 

have .ha. .,.) i^'og/l + ||.). .fVo. /f j,^.. < co. we deduce .ha. «,oo, < oo^ 

Next let a G (— (Tg,0). Lemma O implies that there is a constant c G (0,oo) such that 

^^'^ ~ ^ r , ./^Z . ^^2 (^' + ^'y)''"' dy (132) 

J2 {yP/^log{y)) 

as z — )■ 00. As /2°° (\o \ ))'2 dy < CO, this implies that R{oo) < 00. Finally assume that a < — (Jg. 
Theorem V.54.5 in [32] implies that 

P^ (Zt G dy) ^^ c ^^ exp ( /'' ^^ dz) dy (133) 
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for every z € (0, oo) if there exists a normalizing constant c € (0, oo) such that the right-hand 
side is a probabihty distribution. Due to (jl30p we need to show that 



{mf K + ^ey) "^ = (^^f^) ^ (^' + ^"y) "^ (^^^) 



is integrable over y G (0,oo). This function is bounded over (0, 1] and is of order 0(y °?) as 
y — 7- oo. As a < — Cg, there exists a normahzing constant c such that the right-hand side of ()133p 
is a probabihty distribution. D 

7 Family decomposition of BDREs with immigration 

Let a, G R, (Tf, € (0, oo) and (Tg € [0, oo). In this section we consider the BDRE with immigra- 
tion/emigration which is the solution of the SDEs 

dZt = edt + \alZt dt + ZtdSt + Ja^ZtdW^^^'^ 



re-t^^^^t^^t^yu,^,.r., ^^^^^ 



dSt = adt+ JoUwl''^ 



5t = « ut -|- Y o-a\\ I 

for f > 0. The family decomposition of the BDRE with immigration will be a corollary of the 
family decomposition of Feller's branching diffusion with immigration. For this, we first need to 
generalize Proposition [2] to include immigration. 

Lemma 16. Assume a,6 £ IR., ab £ (0,oo) and a^ G [0, oo). Let {Ft)t>Q be a weak solution of 

dFt = —dt + ^/FtdW^^''^ (136) 

Ae) 



for t S [0,oo) and let St := at + aeW^ for t G [0, oo) be independent of {Ft)t>o- Moreover, 
define the random time (T(t))j>o through r(i) := L e~ "a'^ds for t G [0, oo). Then 

(F,(,)e^%5<)^^Q (137) 

is a weak solution of ()135|) . 

Proof Fix a,6 £ K, ab e (0,oo) and ae G [0,oo). Define Zt := ^^(4)6"^* for t G [0, 00). Ito's 
lemma together with independence of {Ft)t>o and of {St)t>o imply that 

dZt = dF,^t)e^' 

= e^*dF,(t) + F,^t)e^'dSt + \F,^t)e^'al dt 



e^*^dr(t) + e^y F,(i)dVr;(^) + ZtdSt + \ztal dt 

=^*-^r'(t) dt + ^[^p~Jt^^=dwll^ + ^*^^* + ^^*^' ' 
CT,, la'^e~^^ 

1 o„ . „ .„ nr:: i_ 



(138) 



• dt + -aiZt dt + ZtdSt + ^ aiZt-^==dW^l^^ 
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for t G [0,00). As {T{t))t>o and (1^^ )t>o are independent, the process {Wt)t>o defined through 



Wi 



t •- 



^^ VA^)' 



.w 



dW , ■., t £ [0,00), is a continuous martingale and a Markov process satisfying 



r(s) 



E[W^ 



E 



dW 



(b) 



VTiT) -(«) 



E 



Tjr^dTis) 



r'{s) 



(139) 



for all t G [0, 00). Thus {Wt)t>o is a standard Brownian motion according to Levy's characteriza- 
tion (e.g. Theorem IV. 33.1 of [S^). Moreover {Wt)t>o and (VFj )t>o are independent. Therefore 
P38p implies that {Zt,St)t>o is a weak solution of (jl35|) . D 

Let (Tfo E (0,00) and let {Ft)t>o be the solution of the SDE 

dFt = ^piFtdWt (140) 

for t G [0,00). The associated excursion measure Qp on U is uniquely determined through 



g{x) Qridx) 



lim-E^ 
£-5>oe 



gm)t>o) 



(141) 



for all bounded, continuous functions g: C ([0, 00), [0, 00)) — ?• R, depending on a finite time 
interval and such that there is an e > such that g{x) = for all x G C ([0,oo), [0,oo)) with 
supf>o Xt < £• Such a measure Qp exists according to Theorem 1 in plj . 

The following family decomposition of Feller's branching diffusion is a special case of the 
family decomposition of the Dawson- Watanabe superprocess with immigration (see [32] and |llj). 
Recall the excursion space U from (j44p . 

Lemma 17. Let ^ G [0, 00), a G R and cjfe G (0,oo). Let "Pq ^e ^ Poisson point process on 
[0, 00) X [/ TOt/i intensity measure dy x Q^? and let Ve he an independent Poisson point process 
on [0,00) X U with intensity measure 6dt x Qp. Then the process (-Fj)j>o defined through Fq = x 
and 

Ft-= Yl h<-Xt+ Yl ^*-- (142) 

(j/,x)6Po (s,x)6Pe 

for t G (0, 00) is a weak solution of the SDE 



Fn 



(143) 



dFt = edt + ^alFtdWt, 
for t G [0, 00) and for each x G [0, 00). 

Proof of Theorem\3i Fix cri,,ae G (0, 00) and a G (— 00,— dg]. Define a process (Fj)i>o through 
-Fq = 2; and through 

Ft-= y iy<xxt+ y xt-s (144) 



(s/:X)GP (s,x)eP 



for t G (0,00). Then Lemma [T71 shows that {Ft)t>o is a weak solution of 

dFt = dt+ VF\dWt 



(145) 



for t G [0,00). Lemma [TBI implies that {Zt,St)t>o = {Ff{t)^^^ : St)t>o is a weak solution of the 
SDKs ([15|) . Also {Zt,St)t>o is a solution of (|1^ due to Theorem [71 As the solution of (03]) is 
unique in law, we conclude that {Zt,St)t>o and {Zt,St)t>o have the same distribution. D 
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